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O | Abstract 

The spectral measure plays a key role in the statistical modeling of multivariate extremes. 
I Estimation of the spectral measure is a complex issue, given the need to obey a certain moment 

condition. We propose a Euclidean likelihood-based estimator for the spectral measure which 
P— ] is simple and explicitly defined, with its expression being free of Lagrange multipliers. Our 

estimator is shown to have the same limit distribution as the maximum empirical likelihood 
estimator of J. H. J. Einmahl and J. Segers, Annals of Statistics 37(5B), 2953-2989 (2009). 
r/} Numerical experiments suggest an overall good performance and identical behavior to the 

maximum empirical likelihood estimator. We illustrate the method in an extreme temperature 
^ data analysis. 
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Statistics of extremes. 



1 Introduction 



When modeling dependence for bivariate extremes, only an infinite-dimensional object is flexi- 
ble enough to capture the 'spectrum' of all possible types of dependence. One of such infinite- 
dimensional objects is the spectral measure, describing the limit distribution of the relative size 
of the two components in a vector, normalized in a certain way, given that at least one of them 
is large; see, for instance, Kotz and Nadarajah ( 2000[ §3) and Beirlant et al. (2004 §8-9). The 



normalization of the components induces a moment constraint on the spectral measure, making 
its estimation a nontrivial task. 
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In the literature, a wide range of approaches has been proposed. Kotz and Nadarajah (2000 



§2-3) survey many parametric models for the spectral measure, and new models continue to be 



invented ( Cooley et al. 2010 Boldi and Davison 2007 Ballani and Schlather 2011). Here we are 



mostly concerned with semiparametric and nonparametric approaches. Einmahl and Segers (2009) 



propose an enhancement of the empirical spectral measure in Einmahl et al. (2001) by enforcing 



the moment constraints with empirical likelihood methods. A nonparametric Bayesian method 



based on the censored-likclihood approach in Lcdford and Tawn ( 1996 1 is proposed in Guillottc 



et al. (2011) 



In this paper we introduce a Euclidean likelihood-based estimator related with the maximum 



empirical likelihood estimator of |Einmahl and Segers ( 2009 1 . Our estimator replaces the empirical 
likelihood objective function by the Euclidean distance between the barycenter of the unit simplex 



and the vector of probability masses of the spectral measure at the observed pseudo-angles ( Owen 



1991 2001 Crepet et al. 20091. This construction allows us to obtain an empirical likelihood-based 



estimator which is simple and explicitly defined. Its expression is free of Lagrange multipliers, 
which not only simplifies computations but also leads to a more manageable asymptotic theory. We 
show that the limit distribution of the empirical process associated with the maximum Euclidean 
likelihood estimator measure is the same as the one of the maximum empirical likelihood estimator 
in 



Einmahl and Segers (2009). Note that standard large-sample results for Euclidean likelihood 



methods (Xu 1995 Lin and Zhang 2001 ) cannot be applied in the context of bivariate extremes. 



The paper is organized as follows. In the next section we discuss the probabilistic and geometric 
frameworks supporting models for bivariate extremes. In Section[3]we introduce the maximum Eu- 
clidean likelihood estimator for the spectral measure. Large-sample theory is provided in Section^ 
Numerical experiments are reported in Section [5] and an illustration with extreme temperature 
data is given in Section [6] Proofs and some details on a smoothing procedure using Beta kernels 
are given in the Appendix |A"| and |B| respectively. 



2 Background 

Let (Xi, Yi), (X 2 , Y 2 ), . . . be independent and identically distributed bivariate random vectors with 
continuous marginal distributions Fx and Fy. For the purposes of studying extremal dependence, 
it is convenient to standardize the margins to the unit Pareto distribution via X* = 1 /{I— Fx (Xi)} 
and Y* = 1/{1 — Fy(Yi)}. Observe that X* exceeds a threshold t > 1 if and only if X t exceeds 
its tail quantile F^ x (l — 1/i); similarly for Y*. The transformation to unit Pareto distribution 
serves to measure the magnitudes of the two components according to a common scale which is 
free from the actual marginal distributions. 



Pickands' representation theorem ( Pickandsl 19811 asserts that if the vector of rescaled, com 
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poncntwise maxima 

M* = ( — max X* , - max Y* ) , 
yn *=i,...,n n i=i...,n y 

converges in distribution to a non-degenerate limit, then the limiting distribution is a bivariate 
extreme value distribution G with unit-Frechet margins given by 

G(x,y) = exp j - 2 J max - — — ^ dif(w)|, x,y > 0. (1) 

The spectral (probability) measure H is a probability distribution on [0, 1] that is arbitrary apart 
from the moment constraint 

/ wdH(w) = 1/2, (2) 
J [OA] 

induced by the marginal distributions G(z,oo) — G(oo, z) = exp(— 1/z) for z > 0. 

The spectral measure H can be interpreted as the limit distribution of the pseudo-angle Wi = 
X* I (X* + Y*) given that the pseudo-radius Ri = X* + Y* is large. Specifically, weak convergence 
of M* to G is equivalent to 

Pr[Wi £ ■ | Ri >t] ^H(-), t^oo. (3) 

The pseudo-angle Wj is close to or to 1 if one of the components X* or Y* dominates the other 
one, given that at least one of them is large. Conversely, the pseudo-angle Wi will be close to 
1/2 if both components X* and Y* are of the same order of magnitude. In case of asymptotic 
dependence, G(x,y) = exp(— l/x — 1/y), the spectral measure puts mass 1/2 at the atoms and 
1, whereas in case of complete asymptotic dependence G{x,y) = exp{— l/min(a;, y)}, the spectral 
measure H reduces to a unit point mass at 1/2. 

Given a sample {X\, Yi), . . . , (X n , Y n ), we may construct proxies for the unobservable pseudo- 
angles Wi by setting 

W = X*/(X* + Y*), Ri = X* + Y* , 

where X* = 1/{1 - F X {X % )} and Y* = 1/{1 - F Y (Yi)} and where F x = F x . n and F Y = A> are 
estimators of the marginal distribution functions Fx and Fy ■ A robust choice for Fx and Fy is 
the pair of univariate empirical distribution functions, normalized by n + 1 rather than by n to 
avoid division by zero. In this case, X* and Y* are functions of the ranks. 

For a high enough threshold t = t n , the collection of angles {Wi : i £ K} with K — K n = {i = 
1, . . . , n : R\ > t} can be regarded as if it were a sample from the spectral measure H. Parametric 
or nonparamctric inference on H may then be based upon the sample {Wi : i £ K}. Nevertheless, 
two complications occur: 

1 . The choice of the threshold t comes into play both through the rate of convergence in (|3| 
and through the effective size \K\ = k of the sample of pseudo-angles. 
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2. The standardization via the estimated margins induces dependence between the pseudo- 
angles, even when the original random vectors (X,, Yj) are independent. 

For the construction of estimators of the spectral measure, we may thus pretend that {Wi : i G 
1} constitutes a sample from H. However, for the theoretical analysis of the resulting estimators, 
the two issues above must be taken into consideration. Failure to do so will lead to a wrong 
assessment of both the bias and the standard errors of estimators of extremal dependence. 



3 Maximum Euclidean Likelihood Estimator 



We propose to use Euclidean likelihood methods (Owen 2001 pp. 63-66) to estimate the spectral 



measure. Let w%, . . . , Wk € [0, 1] be a sample of pseudo-angles, for example the observed values of 
the random variables Wi, i G K, in the previous section, with k — \K\. The Euclidean loglikelihood 
£e ratio for a candidate spectral measure H supported on {wi, . . . , Wk} and assigning probability 
mass Pi = H({wi}) to Wi is formally defined as 

The Euclidean likelihood ratio can be viewed as a Euclidean measure of the distance of p = 
(pi, . . . ,pk) to the barycenter (& ,...,& ) of the (k — l)-dimensional unit simplex. In this 
sense, the Euclidean likelihood ratio is similar to the empirical loglikelihood ratio 

k 

*(p) = X>g(*Pi), 

1=1 

which can be understood as another measure of the distance from p to . . . , fc" 1 ). Note that 

£e(p) results from £(p) by truncation of the Taylor expension log(l + x) — x — x 2 /2 + • • • and the 
fact that pi + ■ ■ ■ + pk = 1, making the linear term in the expansion disappear. 

We seek to maximize £e(p) subject to the empirical version of the moment constraint Q. Our 
estimator H for the distribution function of the spectral measure is defined as 

k 

H(w)=Y,Pir(w i <w) ) we [0,1], (4) 

i=l 

the vector of probability masses p = {pi, . . . ,pk) solving the optimization problem 

max -iEjUte*- 1 ) 2 
s-t- EtiPi = l ( 5 ) 

ELi w iPi = V 2 - 

This quadratic program with linear constraints can be solved explicitly with the method of La- 
grange multipliers, yielding 

Pi= l{l~{w-l/2)S- 2 {w t ~w)}, i=l,...,k, (6) 
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where w and S 2 denote the sample mean and sample variance of W\ , . . . , Wk, that is, 



w 



1 

EE"*. 



1 fe 



»=1 i=l 

The weights could be negative, but our numerical experiments suggest that this is not as 



problematic as it may seem at first sight, in agreement with Antoine et al. (2007 1 and Crepet 



et al. (2009 1, who claim that the weights pi are nonnegative with probability tending to one. 



The second equality constraint in ^ implies that H satisfies the moment constraint Q, as 
Jj Q j, wdH(w) — Wipi — 1/2, which can be easily verified directly. 

The empirical spectral measure estimator of Einmahl et al. (2001 1 and the maximum empirical 
likelihood estimator of Einmahl and Segers ( 2009 ) can be constructed as in |5| through suitable 
changes in the objective function and the constraints. The empirical spectral measure H(w) = 
^2 i Pi I(wi < w) solves the optimization problem 

max Y,i=i l °gPi 

yielding pi — l/k for i = 1, ...,k. In contrast, the maximum empirical likelihood estimator 
H(w) = ^2 i pi I(wi < w) has probability masses given by the solution of 



PGR* 



(7) 



. A', 



J2i=l w iPi = V 2 - 
By the method of Lagrange multipliers, the solution is given by 

..1 1 . _ 

Pi ~ k 1 + X(w t - 1/2) ' 

where A € K is the Lagrange multiplier associated to the second equality constraint in Q , defined 
implicitly as the solution to the equation 

if Wj-1/2 



^ 1 + X( Wi - 1/2) 



see also Qin and Lawless ( 1994 ) . 



4 Large- Sample Theory 

The maximum Euclidean likelihood estimator H in Q can be expressed in terms of the empirical 
spectral measure H given by 

1 k 

H(w) = -^2l(Wi<w), we [0,1]. 
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Indeed, w and S 2 are just the mean and the variance of H, and the expression ([6| for the weights 
Pi can be written as 

**l{v) = \-{w -1/2) S-\v-w), «e[0,l]. 
dH 

Integrating out this 'likelihood ratio' over v £ [0,w] yields the identity 

#(«>) = ($(H))(tu), we [0,1], 

where the transformation <!> is defined as follows. Let D$ be the set of cumulative distribution 
functions of non-degenerate probability measures on [0,1]. For F £ D$, the function &(F) on 
[0, 1] is defined by 

(<!>(F))(w) = F(w)-(hf-1/2)<t f 2 f (v-fi F )dF(v), we [0,1]. 

J[0,w] 

Here jip — Jj Q ^ vdF(v) and <jp — Jjoi]^ ~ ^f) 2 dF(v) denote the mean and the (non-zero) 
variance of F. 

We view D$ as a subset of the Banach space ^°°([0, 1]) of bounded, real-valued functions on 
[0, 1] equipped with the supremum norm || • W^. The map $ takes values in £°°([0, 1]) as well. 



Weak convergence in ^°°([0, 1]) is denoted by the arrow 'W and is to be understood as in van der 
Vaart and Wellner| ( p96| . 



Asymptotic properties of the empirical spectral measure together with smoothness properties 
of $ lead to asymptotic properties of the maximum Euclidean likelihood estimator: 

• Continuity of the map $ together with consistency of the empirical spectral measure yields 
consistency of the maximum Euclidean likelihood estimator (continuous mapping theorem) . 

• Hadamard differentiability of the map <I> together with asymptotic normality of the empirical 
spectral measure yields asymptotic normality of the maximum Euclidean likelihood estimator 
(functional delta method). 

The following theorems are formulated in terms of maps H n taking values in D$ . The case to 
have in mind is the empirical spectral measure H n (w) — fc„ X^gif I(W% ^ w ) with {Wi : i £ K n } 
and k n = \K n \ as in Section|2] In Theorem 3.1 and equation (7.1) of Einmahl and Segers (20091, 



asymptotic normality of H n is established under certain smoothness conditions on H and growth 
conditions on the threshold sequence t n . 

Theorem 1 (Consistency). If H n are maps taking values in D$ and if \\H n — -ff||oo — > in outer 
probability for some nondegenerate spectral measure H , then, writing H n = <&(./?„), we also have 
\\H n — H | loo — > in outer probability. 



G 



The proof of this and the next theorem is given in Appendix [A] In the next theorem, the 
rate sequence r n is to be thought of as y/k n . Let C([0, 1]) be the space of continuous, real- valued 
functions on [0,1]. 

Theorem 2 (Asymptotic normality). Let H n and H be as in Theorem^ If H is continuous and 
if 

r n (H n - H) ~» j3, ti V oo , 
in 1]), with < r„ — > oo and with f3 a random element o/C([0, 1]), then also 



r n (H n - H) 



7, 



n — > oo 



with 



7 (w) = P(w) - a H 2 / P(v)dv / (l/2-u)dF(u) 



e [o,i]. 



(8) 



Comparing the expression for 7 in d8j) with the one for 7 in (4.7) in Einmahl and Segers 



(2009), we see that the link between the processes /3 and 7 here is the same as the one between 



the processes /3 and 7 in Einmahl and Segers ( 2009 1 . It follows that tuning the empirical spectral 
measure via either maximum empirical likelihood or maximum Euclidean likelihood makes no 
difference asymptotically. The numerical experiments below confirm this asymptotic equivalence. 



To facilitate comparisons with Einmahl and Segers (2009), note that our pseudo-angle w G [0, 1] 
is related to their radial angle 8 £ [0, w/2] via 9 = arctan{u;/(l — w)}, and that the function / in 
(4.2) in |Einmahl and Segersj ( |2009[ ) reduces to f{9) = (sin - cos 0)/(sin 6 + cos 6) =2w-l. 

How does the additional term er^ 2 (3(v)dv f™(l/2 — v)dH(v) influence the asymptotic 
distribution of the maximum Euclidean/empirical estimator? Given the complicated nature of 



the covariancc function of the process (3, see (3.7) and (4.7) in Einmahl and Segers (20091, it 
is virtually impossible to draw any conclusions theoretically. However, Monte Carlo simulations 
in 



Einmahl and Segers (2009 §5.1) confirm that the maximum empirical likelihood estimator is 



typically more efficient than the ordinary empirical spectral measure. These findings arc confirmed 
in the next section. 



5 Monte Carlo Simulations 

In this section, the maximum Euclidean likelihood estimator is compared with the empirical spec- 
tral measure and the maximum empirical likelihood estimator by means of Monte Carlo simula- 
tions. The comparisons are made on the basis of the mean integrated squared error, 



MISE(-) = E 



{■ -H(w)} 2 dw 



The bivariate extreme value distribution G with logistic dependence structure is defined by 



G(x,y) = exp{-(x 



y 



— l/oc\(x 



n, 



x, y > 0, 
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Figure 1: Examples of trajectories of the empirical spectral measure (solid), the maximum empirical 
likelihood estimator (dashed), and the maximum Euclidean likelihood estimator (dotted). The solid 
grey line corresponds to the true spectral measure H a coming from the bivariate logistic extreme value 
distribution with parameters a, = 0.8 (left) and q„ = 0.4 (right). 



in terms of a parameter a S (0,1]. Smaller values of a yield stronger dependence and the limiting 
cases a — > and a = 1 correspond to complete dependence and exact independence, respectively. 
For < a < 1, the spectral measure H a is absolutely continuous with density 

cLH fM = h a -i _ i)wi - w )y 1 - 1 / a {w- 1 / a + (i - w )-^ a } a ~\ o < w < i. 

aw 2 ^ ' 

Here we consider a Y — 0.8 and a a — 0.4, with stronger extremal dependence corresponding to case 
II. For each of these two models, 1000 Monte Carlo samples of size 1000 were generated. The 
thresholds were set at the empirical quantiles of the radius R given by t = 75%, 75.5%, . . . , 99.5% 
for case I and t — 50%, 50.5%, . . . , 99.5% for case II. The margins were estimated parametrically, 
by fitting univariate extreme value distributions using maximum likelihood. 

In Figure [T] a typical trajectory of the estimators is shown, illustrating the closeness of the 
maximum Euclidean and empirical likelihood estimators. The good performance of the maximum 
Euclidean/empirical spectral measure is confirmed by Figure [2j For larger k (lower threshold t), 
the bias coming from the approximation error in ([3| is clearly visible. 

Numerical experiments in Einmahl and Segers| ( |2009| §5.2) show that the presence of atoms at 
the endpoints and 1 has an adverse effect on maximum Euclidean/empirical likelihood estimates, 
and this finding is further confirmed by Guillotte et al.| (2011 §7.1). Indeed, by construction, the 
pseudo-angles Wi will never be exactly or 1. The empirical spectral measure therefore does not 
assign any mass at and 1, and this situation cannot be remedied by the maximum empirical or 
Euclidean likelihood estimators, having the same support as the empirical spectral measure. 




Figure 2: Logarithms of the mean integrated squared errors of the spectral measure estimates based on 
1000 samples of size 1000 from the bivariate logistic extreme value distribution with parameters ai = 0.8 
(left) and an = 0.4 (right). The solid, dashed, and dotted lines correspond to the empirical spectral 
measure, the maximum empirical likelihood estimator, and the maximum Euclidean likelihood estimator, 
respectively. 



The weights pi of the maximum Euclidean likelihood estimator can be negative. However, 
as can be seen from Figure [3j the weights tend to be positive overall, except for extremely high 
thresholds, with the proportion of negative weights being smaller in case I. This suggests that the 
closer we get to exact independence, the lower the proportion of negative weights. 



50 100 150 200 250 
k 



100 200 300 400 500 
k 



Figure 3: Proportion of negative weights based on 1000 samples of size 1000 from the bivariate logistic 
extreme value distribution with parameters a r = 0.8 (left) and a u = 0.4 (right). 
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6 Extreme Temperature Data Analysis 



6.1 Data Description and Preliminary Considerations 

The data were gathered from the Long-term Forest Ecosystem Research database, which is main- 
tained by LWF (Langfristige Waldokosystem-Forschung) , and consist of daily average meteoro- 
logical measurements made in Beatenberg's forest in the canton of Bern, Switzerland. More 
information on these data can be found at 

\protect\vrule widthOpt\protect\href {http : //www. wsl . chMhttp : //www. wsl . ch} 



and for an extensive study see Ferrez et al. (2011 1. Two time series of air temperature data are 



available: One in the open field and the other in a nearby site under the forest cover. Our aim 
is to understand how the extremes in the open relate with those under the canopy; comparison 
of open-site and below-canopy climatic conditions is a subject of considerable interest in Forestry 



and Meteorology (Renaud and Rebetez 2009 Ferrez et al. 2011 Renaud et al. 20111. The raw 



data are plotted in Figure |4j but before we are able to measure extremal dependence of open air 
and forest cover temperatures we first need to preprocess the data. The preprocessing step is the 



same as in Ferrez et al. (2011 §3.1) and further details can be found in there. 




1998 2000 2002 2004 2006 1998 2000 2002 2004 2006 

Time (in years) Time (in years) 



Figure 4: Daily average air temperatures from the meteorological station in Beatenberg's forest: under 
the forest cover (left) and in the open field (right). 



We consider daily maxima of the residual series that result from removal of the annual cycle in 
both location and scale, and we then take the residuals at their 98% quantile; hence the threshold 
boundary is defined as U = {(x,y) G [0,oo) 2 : x + y = F'^ 1 (0.98) = 105.83}, so that there are 
\U\ = k = 57 exceedances. 
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Figure 5: Scatterplot of air temperature data after transformation to unit Frechet scale; the solid line 
corresponds to the boundary threshold in the log-log scale, with both axes being logarithmic. 
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Figure 6: Empirical estimates and 95% pointwise confidence intervals for x( u ) (left) and x( u ) (right) as 
a function of u £ (0, 1). 



The dependence between open field and forest cover temperatures can be observed in Figure [5j 
where we plot a log- log scale scatterplot of the unit Frechet data, and where we note that after 
log transformation the linearity of the threshold boundary U is perturbed. 
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Spectral measures are only appropriate for modeling asymptotically dependent data. Although 



this issue has been already addressed in Ferrez et al. (2011 1, for exploratory purposes we present 



in Figure [6] the empirical estimates of the dependence coefficients x( u ) an d x( u )> f° r < u < 1, 
defined in |Coles et al. ( 1999[ ) as 

Pr[F x (X) < u,F Y (Y) < u] 21og(l-«) 



X (u) = 2 - log 



logFx[F x (X) < u] ' AW \ogVx[F x {X) >u,F Y (Y) > u] 
Although these plots fail to have a clear-cut interpretation given the large uncertainty entailed in 
the estimation, the point estimates seem to be consistent with asymptotic dependence as already 



noticed by Ferrez et al. (2011 1. 



6.2 Extremal Dependence of Open Air and Forest Cover Temperatures 

We now apply the maximum Euclidean likelihood estimator to measure extremal dependence of 
open air and forest cover temperatures. The estimated spectral measure is shown in Figure [7] All 
weights are positive, i.e. pi > 0, for i = 1, . . . , 57. 

By construction, the estimate of the spectral measure is discrete. A smooth version which still 
obeys the moment constraint ^ can easily obtained by smoothing the maximum Euclidean or 
empirical likelihood estimator with a Beta kernel. Related ideas are already explored in |Hall and| 



Presnell 


(1999 


1 and 


Chen 


(1997 



A cross- validatory procedure was used to select the bandwidth, yielding a concentration pa- 



rameter of v ss 163. Numerical experiments in Warchol (2012) suggest that convoluting empiri- 



cal likelihood-based estimators with a Beta kernel yields a further reduction in mean integrated 



squared error. The Beta kernel even outperforms Chen's kernel (Chen 1999), which is asymp- 



totically optimal under some conditions (Bouezmarni and Rolin 20031, but which is unable to 
conserve the moment constraint. 

From the smoothed spectral measure, we obtain an estimate of the spectral density and plug-in 
estimators for the Pickands dependence function A(w) = 1 — w + 2 J Q H (v) dv, w G [0,1], and the 
bivariate extreme value distribution in (JT|). The estimated spectral density is compared with the 
fit obtained from the asymmetric logistic model 



H, 



i + 1>i + ^2 - {Vh 1/a (i - w) 1 ^" 1 - V> 1/Q_1 } 

{V>l /a (l - H 1/Q + 4 /a w 1/a } 1/1/a - 



e [0,1], 



with parameter estimates -01 = 0.78 (standard error 0.03), V>2 = 0.90 (0.03) and a = 0.30 (0.02). 



The asymmetric logistic model was considered by Ferrez et al. (2011) as the parametric model 
that achieved the "best overall fit." 

In Figure [7] we also plot the smooth spectral measure and corresponding spectral density 
which are obtained by suitably convoluting the empirical Euclidean spectral measure with a Beta 
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Figure 7: Estimates of the spectral measure (left) and the spectral density (right). Left: the dotted line 
corresponds to the empirical Euclidean spectral measure, and the solid line corresponds to its smooth 



version constructed using (10 1. Right: The solid line corresponds to the smooth spectral density obtained 



from the Euclidean likelihood weights using (J9J, and the gray line represents the fit from the asymmetric 
logistic model with (V>i, $ 2 ) = (0.78, 0.90) and p = 0.30. 



kernel as described in (|9| and (10). Since more mass concentrated over 1/2 corresponds to more 
extremal dependence, and more mass concentrated on and 1 corresponds more independence in 
the extremes, a rough interpretation for our context is as follows: The lower the shelter ability 
of the forest, the more mass should be concentrated around 1/2, whereas higher shelter ability 
corresponds to the case where the spectral measure gets more mass concentrated at 1; more mass 
concentrated at suggests relatively more extreme events under the forest cover, suggesting that 
the forest has the ability to retain heat during extreme events. 

In Figure [8] 1 we plot the corresponding Pickands dependence function. More extremal depen- 
dence corresponds to lower Pickands dependence functions, and the deeper these are on the right 
the less frequent are the extreme events under the forest cover relatively to the open field. Our 
analysis suggests that extreme high temperatures under the forest cover are more frequent than 
expected from a corresponding parametric analysis. This somewhat surprising finding is already 



predicted in Ferrez et al. (2011 Fig. 4). The phenomenon may be due to the ability of some forests 
to retain heat, acting like a greenhouse, or it may be connected with the way that other features 



of the forest's structure can alter its microclimate (Renaud et al. 2011 1. Along with the Pickands 
dependence function, we also plot in Figure [HJl the pseudo-angles which provide further evidence 
of a marked right skewness. 

The joint behavior of temperatures in the open and under forest cover can also be examined 
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Figure 8: Estimates of the Pickands dependence function and contours of the bivariate extreme value 
distribution. Left: The solid line represents the smooth Pickands dependence function obtained from the 



Euclidean likelihood weights using (12 1, and the gray line represents the fit from the asymmetric logistic 



model with (^1,^2) = (0.78,0.90) and p = 0.30. The pseudo-angles are presented at the top. Right: The 
solid gray line represents the smooth bivariate extreme value distribution obtained from the Euclidean 
likelihood weights using Jl3p . Air temperature data are plotted on the unit Frechet scale. 



from the estimated bivariate extreme value distribution function plotted in Figure [8]2, which 
was constructed by convoluting the empirical Euclidean spectral measure with a Beta kernel as 



described in (13) 



7 Discussion 

In this paper we propose a simple empirical likelihood-based estimator for the spectral measure, 
whose asymptotic efficiency is comparable to the empirical likelihood spectral measure of Einmahl 



and Segers ( 2009 1 . The fact that our estimator has the same limit distribution as the empirical 



likelihood spectral measure, suggests that a more general result may hold for other members 
of the Crcssic Read class, of which these estimators are particular cases, similarly to what was 



established by Baggerly (19981 in a context different from ours. We focus on the spectral measure 



defined over the Li-norm, but only for a matter of simplicity, and there is no problem in defining 
our estimator for the spectral measure defined over the L p -norm, with p € [1, 00]. For real data 
applications smooth versions of empirical the estimator may be preferred, but these can be readily 
constructed by suitably convoluting the weights of our empirical likelihood-based method with a 
kernel on the simplex. 
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A Proofs 

Proof of Theorem^ Let F n , F e D$. By Fubini's theorem, 



vdF(v) =wF(w) - / F(v)dv, we [0,1], 

[0,w] JO 



[ v 2 dF(v) = l- [ 2vF(v)dv, 
J\o.i] Jo 



'[0 

and similarly for F n . It follows that \\F n — i^Hoo ~~ t implies that \ip n — > [ip, o~ F —> <j f , and 
J [Qw] vdF n (v) -> J [Qw] vdF(v) uniformly in w £ [0, 1]. Hence \\${F n ) - *(i ? )|| 00 -> 0. Therefore, 
the map $ : D$ — > i°°([0, 1]) is continuous. The lemma now follows from the fact that $>(H) = H 
together with the continuous mapping theorem, see Theorem 1.9.5 in |van der Vaart and Wellner| 



Q1996D - □ 
Proof of Theorem^ Write 

fin = r„(H n - H), "f n = r n (H n - H). 

Let D n denote the set of functions / £ £°°([0, 1]) such that H + r^ 1 f belongs to D$. Since 
H + r~ P n = H n takes values in D$, it follows that (3 n takes values in D„. Define g n : D„ — > 
£°°([0,1]) by 

9n (f) = r n mH + r- 1 f)-H}. 

Observe that 

<?„(/?„) = r n {$(H n ) -H} = 7n . 
Further, define the map g : C([0, 1]) -> C([0, 1]) by 

p 1 pW 

(9(/))H = fH-a H 2 f(v)dv (l/2-v)dH(v), we [0,1]. 
Jo Jo 



15 



A straightforward computation shows that if /„ e D n is such that ||/ n — /||<x> ~^ for some 
/ € C([0, 1]), then \\g n (f n ) — <7(/)||oo — * 0; note in particular that /„(!) = and thus also /(l) = 0. 



The extended continuous mapping theorem (van der Vaart and Wellner 1996 Theorem 1.11.1) 
implies that 

In = 9n(/3 n ) -> g{0) = 7, n -> OO, 

as required. Note that we have actually shown that $ is Hadamard differentiable at iZ tangentially 
to C([0, 1]) with derivative given by <I>^ = g. The result then also follows from the functional delta 
method. □ 



B Beta-Kernel Smoothing of Discrete Spectral Measures 

We only consider the case of the empirical Euclidean spectral measure using a Beta kernel, but 
the same applies to the empirical likelihood spectral measure by replacing the pi with pi in ([3]). 
The smooth Euclidean spectral density is thus defined as 

k 

h(w) = } j p i l3{w;WiL',(l-Wi)v}, we (0,1), (9) 

i=l 

where v > is the concentration parameter (inverse of the squared bandwidth, to be chosen 
via cross-validation) and /3(w;p,q) denotes the Beta density with parameters p, q > 0. The 
corresponding smoothed spectral measure is defined as 

r w ~ 

H(w)= h{v)dv = y2piB{w; Wi v,{l-Wi)v}, we[0,l], (10) 

where B(w;p, q) is the regularized incomplete beta function, with p, q > 0. Since 

/ wh(w)dw = y2pi\ ■ — - q I = = 1/2, 

Jo L VWi + 1/(1 - Wi) J ^ 



(11) 



the moment constraint is satisfied. Plug-in estimators for the Pickands dependence function and 
the bivariate extreme value distribution follow directly from 

A(w) = 1 -w + 2 Vp, \ B{u;w l v,{l-w l )v}du 1 we [0,1], (12) 
i=i Jo 

G(x, y) = exp j-^y^j?i J max 1 - - ^J j3{u;WiU, (1 - it>i)^}du| , a;,y>0. (13) 
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